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Abstract 

We investigate solutions of backward stochastic differential equations (BSDE) 
with time delayed generators driven by Brownian motions and Poisson ran- 
dom measures, that constitute the two components of a Levy process. In this 
new type of equations, the generator can depend on the past values of a solu- 
tion, by feeding them back into the dynamics with a time lag. For such time 
delayed BSDE, we prove existence and uniqueness of solutions provided we 
restrict on a sufficiently small time horizon or the generator possesses a suffi- 
ciently small Lipschitz constant. We study differentiability in the variational 
or Malliavin sense and derive equations that are satisfied by the Malliavin 
gradient processes. On the chosen stochastic basis this addresses smoothness 
both with respect to the continuous part of our Levy process in terms of the 
classical Malliavin derivative for Hilbert space valued random variables, as well 
as with respect to the pure jump component for which it takes the form of an 
increment quotient operator related to the Picard difference operator. 
Keywords: backward stochastic differential equation, time delayed generator, 
Poisson random measure, Malliavin's calculus, canonical Levy space, Picard 
difference operator. 



1 Introduction 

Introduced in [2D], backward stochastic differential equations have been thoroughly 
studied in the literature during the last decade, see [13] or [15] and references therein. 
Viewed from the perspective of Peng who interprets their key structural feature as 
a nonlinear conditional expectation, the close link to the stochastic calculus of vari- 
ations or Malliavin's calculus becomes apparent. In fact, in a Clark-Ocone type 
formula, the control component of the solution pair of a BSDE with a classical glob- 
ally Lipschitz generator without time delay on a Gaussian basis turns out to be the 
Malliavin trace of the other component, see Proposition 5.3 in [13j or Theorem 3.3.1 
in [12]. Not only this observation attributes an important role to Malliavin's calculus 
in the context of stochastic control theory and BSDE. As the simplest example, let 
us recall that hedging strategies in complete market models corresponds to Malliavin 
derivatives of wealth processes, see [Tf]. The fine structure and sensitivity proper- 
ties of solutions of BSDE or systems of forward and backward stochastic differential 
equations have been approached by means of the stochastic calculus of variations 
(see [2] and [Ij), and applied to provide explicit descriptions of delta hedges of insur- 
ance related financial derivatives in [3]. Let us mention that Malliavin's calculus has 
been applied to prove regularity of trajectories and thus to provide a first numerical 
scheme for BSDE with generators of quadratic growth, see for instance \TU\. More 
generally, it has been established as a key tool in the numerics of control theory 
and mathematical finance, for instance to enhance the convergence speed of dis- 
cretization schemes for solutions of BSDE, see jllj . |18| . BSDE have proved to be 
an efficient and powerful tool in a variety of applications in stochastic control and 
mathematical finance. In all of these applications, variational smoothness of their 
solutions is fundamental for describing their properties. 

In this spirit, and with the aim of clarifying smoothness in the sense of the 
stochastic calculus of variations and related properties of BSDE in a more general 
setting, in this paper we study the equations with dynamics given for t G [0, T] by 

Y{t) = i+ j f{s,Ys,Zs,Us)ds- f Z{s)dW{s)- j U{s,z)M{ds,dz). 

An equation of this type will be called BSDE with time delayed generator. It is 
driven by a Levy process, the components of which are given by a Brownian mo- 
tion and a Poisson random measure. In this new type of equations, a generator 
/ at time s depends in some measurable way on the past values of a solution 
{Ys, Zs, Us) = (y{s + u), Z{s + u), U{s + u, .))_T<M<o- Very recently, time delayed 
BSDE driven by Brownian motion and with Lipschitz continuous generators have 
been investigated for the the first time in [8], and in more depth in [10]. We would 
like to refer the interested reader to the accompanying paper [TU], where existence 



and uniqueness questions are treated, and examples given in which multiple solu- 
tions or no solutions at all exist. Further, several solution properties are investigated, 
including the comparison principle, measure solutions, the inheritance property of 
boundedness from terminal condition to solution, as well as the BMO martingale 
property for the control component. We would like to point out that all results from 
|10] can be extended and proved in the setting of this paper. 

Our main findings are the following. First, we prove that a unique solution ex- 
ists, provided that the Lipschitz constant of the generator is sufficiently small, or the 
equation is considered on a sufficiently small time horizon. This is the extension of 
Theorem 2.1 from [1^ to be expected. Secondly, we establish Malliavin's differentia- 
bility of the solution of a time delayed BSDE, both with respect to the continuous 
component of the Levy process, which coincides with the classical Malliavin deriva- 
tive for Hilbert-valued random variables, as well as with respect to the pure jump 
part, in terms of an increment quotient operator related to Picard's difference op- 
erator. We prove that the well-known connection between (Z, U) and the Malliavin 
trace of Y still holds in the case of time delayed generators. 

BSDE without time delays and driven by Poisson random measures have already 
been thoroughly investigated in the literature, see [5], [6] or [23]. But contrary to the 
case with a Gaussian basis, smoothness results in the sense of Malliavin's calculus 
have not been established yet in a systematic way. To the best of our knowledge, 
only in [7j, variational differentiability of a solution of a forward-backward SDK with 
jumps with respect to the Brownian component is considered while differentiability 
with respect to the jump component is neglected. 

We would like to emphasize that backward stochastic differential equations with 
time-delayed generators arise in financial and insurance problems dealing with pric- 
ing, hedging, risk management and optimal control, see the working paper |0. For 
instance, they are encountered in the context of the optimal liquidation problem 
of large trader's positions. A related optimal control problem in terms of BSDE 
exhibits generators in which the delayed feedback of the large trader's actions on 
the price dynamics take the form of a delay effect in the sense considered in this 
paper. As explained at the beginning of this section, Malliavin's calculus plays a 
fundamental role in mathematical finance and optimal control. We believe that the 
results concerning Malliavin's differentiability obtained in this paper are as impor- 
tant to describe parameter sensitivity properties of financial derivatives in this more 
general setting as they are in j3] for generalizing the Black-Scholes delta hedge to 
incomplete markets in a purely probabilistic approach via BSDE. 

This paper is structured as follows. Section 2 deals with the existence and unique- 
ness problem. In Section 3 we survey concepts of the canonical Levy space and 
variational differentiation, and prove some technical lemmas. The main theorem 



concerning Malliavin smoothness of a solution, and the interpretation of the latter 
in terms of a Malliavin trace is proved in Section 4. 

2 Existence and uniqueness of a solution 

We consider a probability space (fi, J-", P) with a filtration F = {J^t)o<t<T, where 
T < oo is a finite time horizon. We assume that the filtration F is the natural 
filtration generated by a Levy process L := {L{t), < t < T) and that J-q contains 
all sets of P-measure zero, so that the usual conditions are fulfilled. As usual, by 
B{X) we denote the Borel sets of a topological space X, while A stands for Lebesgue 
measure. 

It is well-known that a Levy process satisfies the Levy-Ito decomposition 

L{t) =at + aW{t) + / / zN{ds, dz) + / z{N{ds, dz) - u{dz)ds), 

Jo J\z\>l Jo Jo<\z\<l 

for < t < T, with a e R,a > 0. Here W := {W{t),0 < t < T) denotes a 
Brownian motion and A^ a random measure on [0, T] x (R — {0}), so that W and A^ 
are independent. The random measure A^ 

N{t, A) = tl{0 < s < t; AL(s) G A}, < t < T, A G -B(M - {0}), 

counts the number of jumps of a given size. It is called Poisson random measure 
since, for t G [0, T] and a Borel set A such that its closure does not contain zero, 
A^(t, A) is a Poisson distributed random variable. The a-finite measure z/, defined 
on i3(R — {0}), appears in the compensator \® v oi the random measure A^. The 
compensated Poisson random measure (or martingale-valued measure) is denoted 
hj N{t,A) = N{t,A)-tv{A), t G [0,T],Ag i3(R-{0}). In this paper we deal with 
the random measure 

M{t,A) = I I zN{ds,dz) 

Jo J A 

= zN{ds, dz)- / zu{dz)ds, 0<t<T,Ae B{R - {0}). 

Jo J A Jo J A 

It can be considered as a compensated compound Poisson random measure as, for 
a fixed t G [0,T] and a Borel set A the closure of which does not contain zero, 
Jg J^ zN{ds, dz) is a compound Poisson distributed random variable. Finally, we 
introduce the a-finite measure 



m{A) = [ z^u{dz), AeB(R- {0}). 
Ja 



For details concerning Levy processes, Poisson random measures and integration 
with respect to martingale-valued random measures we refer the reader to Chapter 



2 and Chapter 4 of |1]. 

Let us now turn to the main subject of this paper. We study solutions {Y, Z, U) : = 
(F(t), Z(t), f/(t, z))o<i<T,2G(iR-{o}) of a BSDE with time delayed generator, the dy- 
namics of which is given by 

Y{t) = i+ I f{s,Ys,Z,,Us)ds 



t 

T pT 



-f Z{s)dW{s)- f f U{s,z)M{ds,dz), < t < T. (2.1) 

The generator / depends on the past values of the solution, fed back into the system 
with a time delay, denoted by Yg := {Y{s + v))^t<v<o, ^s •= iZ{s + v))-t<v<o and 
Us := (f/(s + v,.))^T<v<o,0 < s < T. We always set Z{t) = U{t,.) = and 
Y{t) = Y{0) for t < 0. Note that the measure M, not A^, is taken to drive the jump 
noise. The reason for this is that we adopt the concepts of Malliavin calculus on the 
canonical Levy space from [21], which is formulated in terms of multiple stochastic 
integrals with respect to M. 

We shall work with the function spaces of the following definition. 

Definition 2.1. 1. Let L^j.(]R) denote the space of measurable functions z : 

[— T, 0] — )■ M satisfying 

1-0 

/ \z{t)\^ dt < oo. 

2. Let L'^_rp^{E) denote the space of product measurable functions u : [— T, 0] x 
(R - {0}) -^ M satisfying 

I- 

2 



|-u(t,2;)| m{dz)dt < oo. 
T Jm-{o} 

3. Let L'^y(M) denote the space of bounded, measurable functions y : [— T, 0] — > M 

such that 

sup \y[t)\ < oo. 
te[-T,o] 

4- Let L^(]R) denote the space of J^T-f^^osurable random variables ^ : fi — )■ M 
which fulfill 

E[|^|^] <oo. 

5. Let EI|.(R) denote the space of predictable processes Z : fi x [0, T] — )■ R such 
that 

t-T 

E[ / \Z{t)\^dt\ <cx). 
Jo 



6. Let E[|,^(]R) denote the space of predictable processes f/ : fi x [0,T] x (M — 
{0}) —7- M satisfying 

-T 



E[ / \U{t, z)\'^m{dz)dt] <oo. 

Jo jR-iO} 



'O JM-jO} 

7. Finally, /etS|.(R) denote the space of¥-adapted, product measurable processes 

F : fi X [0, T] ^ M satisfying 

E[ sup \Y{t)\^] <oo. 

iG[0,T] 

The spaces EI|.(R),E[|,^(M) and S|,(M) are endowed with the norms 

\\Z\\l = E[ r e(^^\Z{t)\'dt], 
^ Jo 

\\U\\l, = E[ [ [ e^'\U{t,z)\^m{dz)dt], 
"^•"^ Jo Jm-io) 



'0 JM-jO} 

\Y\\l, = E[ sup e^*|r(t)|'], 
^ te[o,r] 



with some /3 > 0. 



Predictability of Z means measurability with respect to the predictable cr-algebra, 
which we denote by V, while predictability of U means measurability with respect to 
the product V<S)B{R-{0}). In the sequel let us simply write §2(]R) x lf{R) x H^^(M) 
instead of §^(M) x e^(M) x e^^„(M). 

We start with establishing existence and uniqueness of a solution of (12. ip under 
the following hypotheses: 



(Al) the terminal value ^ G L^( 

(A2) m is a finite measure, i.e. /ig_rQi z'^i'{dz) < oo, 

(A3) the generator / : fi x [0,T] x L^^(R) x L^y(M) x ^^^^^(R) ^ M is prod- 
uct measurable, F-adapted and Lipschitz continuous in the sense that for a 
probability measure a on ([— T, 0], B{[—T, 0])) and with a constant K > 

\f{u, t, yt, Zt, ut) - f{uj, t, yt, Zu ^t) T 
< k( \y{t + v)-y{t + v)\^a{dv)+ I \z{t + v) - z{t + v)\^ a{dv) 

+ / \u{t + V, z) — u{t + V, z)\'^m{dz)a{dv)], 

J-T iR-{0} ^ 

holds for P® A-a.e. (w, t) e i7 x [0, T], for any {yt, Zt, Ut), {yt, Zt, Ut) G L°°^(]R) x 
Ll^{R) X L2 (M), 



(A4) E[J^\f{t,0,0,0)\^dt] <oo, 

(A5) f{u, t, ., .,.) = OioT u eQ,t < 0. 

For convenience, in the notation of / the dependence on u is omitted and we write 
f{t, .,.,.) for f{u, t, ., ., .) etc. We remark that f{t, 0, 0, 0) in (A4) should be under- 
stood as the value of the generator f{t,yt,Zt,Ut) at yt = Zt = Ut = 0. We would 
like to point out that assumption (A5) in fact allows us to take Y{t) = Y{0) and 
Z{t) = U{t, .) = for t < as a solution of (12. ip . Finally, let us recall that under 
(A2) and for an integrand U G E[^(R), the stochastic integral with respect to the 
martingale-valued measure M 



I I U{s,z)M{ds,dz), 0<t<T, 

Jo JR-IOj 



/O JR-{0} 

is well-defined in the Ito sense, see Chapter 4.1 in [4J. 

First let us notice that for {Y, Z, U) E E>'^{R) x H2(]R) x H^(]R) the generator is 
well-defined and integrable as a consequence of 

\fit,Yt,Zt,Ut)\''dt<2 [ \f{t,0,0,0)\^dt + 2K{[ [ \Y{t + v)\''a{dv)dt 

Jo Jo J-T 

pT pO pT nO r 

+ / \Z{t + v)\'^a{dv)dt+ / / / \U{t + V, z)\^m{dz)a{dv)dt) 

Jo J-T Jo J~T iR-{0} 

j-T j-0 j-T+v 

= 2/ \f{t,0,0,0)\'^dt + 2K / \Y{w)\'^dwa{dv) 

Jo J-T Jv 

/O j-T+v 
/ \Z{w)\'^dwa{dv) 

/O pT+v p 
/ / \U{w,z)\'^m{dz)dwa{dv) 

T Jv Jr-{0} 

< 2 /" |/(t,0,0,0)prft + 2K(T sup |y(w;)p 
Jo we[o,T] 

+ \Z{w)\'^dw+ / \U{w, z)\'^m{dz)dw) <oo, ¥ - a.s., (2.2) 

Jo Jo Jr-{o} 

where we apply (A3), Fubini's theorem, use the assumption that Z{t) = U(t, .) = 
and Y{t) = Y{0) for t < and the fact that the measure a is a probability measure. 
The main theorem of this section is an extension of Theorem 2.1 from \1U\ . 
Although the extension is quite natural, the proof is given for completeness and 
convenience of the reader. The key result follows from the following a priori esti- 
mates. 



Lemma 2.1. Let {Y, Z, U), {Y, Z, U) e §2(M) x B.'^(R) x H^^M) denote solutions of 
02.1 p with corresponding parameters {^, f) and {^, f) which satisfy the assumptions 
(Al)-(A5). Then the following inequalities hold 

\\z-z\\l, + \\u-u\\l.^ 

< e^^E[|e-eT]+4E[ / e^V{t,Yt,Zt,Ut)-f{t.YuZt,Ut)?dt], (2.3) 

P Jo 

\\y-y\\l^ (2.4) 

< Se'^^Efle-ll'j+STEf / e^'\f{t,Yt,Zt,Ut)-f{t,Yt,Zt,Ut)\Ht]. 

Jo 

Proof: 

The inequality (12.31) follows by a straightforward extension of Lemma 3.2.1 from 
[15] . by only adding an additional stochastic integral with respect to M. In order to 
prove the second inequality, first notice that for t G [0,T] 

Y{t) - Y{t) = E[e - I + y {f{s, n, Z,, Us) - f{s, n, Z,, Us))ds\:Ft\ , 
and 

e^'\Y{t)-Y{t)\ 

< e2^E[|e-e|l-^t]+E[ / e^-'\fis,Ys,Zs,Us)-fis,Ys,Zs,Us)\ds\Tt], 

Jo 

hold P-a.s.. Doob's martingale inequality and Cauchy-Schwarz' inequality yield the 
second estimate. The reader may also consult Proposition 2.2 in [5.J or Proposition 
3.3 in [6], where similar estimates for BSDE with jumps are derived. D 



Theorem 2.1. Assume that (A1)-(A5) hold. For a sufficiently small time horizon 
T or for a sufficiently small Lipschitz constant K of the generator f , more precisely 
if for some (3 > 

1 /■" 

6{T, K, /3, a) := (8T + -)K / e'^'^aidv) max{l, T} < 1, 

P J -T 

the backward stochastic differential equation (12. ip has a unique solution [Y, Z, U) G 
§2(R) X H2(M) X M^ 



Proof: 

We follow the classical Picard type iteration scheme (see Theorem 2.1 in [13j or 



Theorem 3.2.1 in p!5|) to prove existence and uniqueness of a solution. 

Let Y\t) = Z^{t) = U^{t,z) =0, {t,z) G [0,T] x (R-{0}). 
Step 1) We show that the recursive definition 



rn+l 



- f Z''^\s)dW{s)- f f U''^\s,z)M{ds,dz), 0<t<T, (2.5) 



makes sense. More precisely, we show that given (F",Z",f/") G S^(]R) x EI^(M) x 
H^(R), equation ([23]) has a unique solution (r"+i, Z"+i, f/"+i) G §2(M) x H2(M) x 



^m 



2 

Applying inequality (12.21) . we can conclude that 



E[[ \f{t,Y;\Z:^,UnM < 2E[/ \f{t,0,0,0)\'dt] 
Jo Jo 

+2ir(T||r"||s2 + ||Z"||h2 + ||f/"||H2^) < oo. 

As in the case of BSDE without time delays, the martingale representation, see The- 
orem 13.49 in [14J, provides a unique process Z"'^^ G HI^(M) and a unique predictable 
process U"-^^ satisfying 



E[ /" /" IW'+^t, z)\^u{dz)dt\ <oo, 
Jo Jr-w\ 



so that 



e+/ f{t,Y,\Z^,Ul^)dt = ¥.[i+ I fit,Y,\Z^,Undt] 
Jo Jo 



T pT 



+ Z''+\t)dW{t)+ / U''+\t,z)N{dt,dz), ¥-a.s. 

Jo Jo Jr-{o} 

For (t, z) G [0,T] x (M - {0}) we get t/^+^i, z) = T^^^Shil ^ e^(M), and have the 
required representation 

e+Z" f{t,Y,-,Z^,Undt = E[^+ [ f{t,Y,\Z^,Undt] 
Jo Jo 

+ [ Z''+\t)dW{t)+ f f U''+\t,z){t)M{dt,dz), F-a.s.. 
Jo Jo Jr-{o} 

Finally, we take y"+i as a progressively measurable, cadlag modification of 

Y-+\t)iu) = E[^ + J f{s,Y:,Z:,U:))ds\Tt\, uen,te[0,T]. 



10 



Similarly as in Lemma [2. ![ Doob's martingale inequality, Cauchy-Schwarz' inequal- 
ity and the estimates ([22D yield that y"+i G S'^{R). 

Step 2) We prove the convergence of the sequence ("K", Z", [/"■) in S^(]R) x H^(]R) x 

H^(R). 

The estimates (12. 3p and (12. 4p provide the inequality 

<{8T+^)E[I e^y{t,Y,^,z^,un - f{t,Yr\zr\ur')M. (2.6) 

By applying the Lipschitz condition (A3), Fubini's theorem, changing variables and 
using the assumption Vn > F"(s) = F"(0) and Z'^(s) = U"'{s, .) = for s < 0, we 
can derive 



Jo 

< KE[ [ e^* / \Y''{t + v)-Y''-\t + v)\^aidv)dt 

Jo J-T 

T pO 



e^* / \Z''{t + v)- Z^'^^t + v)Ya{dv)dt 

IQ J-T 

+ f e^' [ [ \U''it + v,z)-U''-\t + v,z)\^m{dz)a{dv)dt\ 

Jo J~T Jr~{0} 

KE[ f e-^" / e^(*+")|r"(t + i;)-r"-i(t + t;)|2rfto(dt;) 

J~T Jo 

J-T Jo 

+ [ e-^" [ I e^^'+''^\U''{t + v,z)-U''-\t + v,z)\''m{dz)dta{dv)\ 

J-T Jo JR-jO} 

/O j-T+v 

e"^" / e^"'|r"(^) - Y''~\w)\''dwa{dv) 
-T Jv 

/O pT+v 

^-fiv / e'^'"|Z"H - Z''-^{w)\^dwa{dv) 
T Jv 

j-T+v j- 

e-^'" / e^'"\U''{w,z)-U''-\w,z)\''m{dz)dwa{dv)\ 

-T Jv Jn-{0} 

j-0 

< K e "^ a{dv)[l\\Y -Y \\^:, + \\Z -Z \y2 + \\U -U ||u2 j- 

J-T '" 

(2.7) 



11 



From (12. 6p and (12. 7p . we obtain 



F ^^ ll§2 + 11^ -^ IIh2 + Ir -'^ IIh2, 

< 6{T,K,P,a){\\Y'^ - Y^-'Wl + ||Z" - Z^-'\\l, + ||t/" - t/"-i||^^),(2.8) 

with 

1 /"^ 
6{T, K, /3, a) = (8T + -)K / e-^^a^dv) max{l, T}. 

For (3 = ^ we have 

S{T,K,/3,a) < 9TKemax{l,T}. 

For sufficiently small T or sufficiently small K, the inequality (12. 8p provides a unique 
limit (Y, Z, U) G §2(m) xe2(R) xH^(M) of the converging sequence (F", Z", f/")„6N, 
which satisfies the fixed point equation 

Y{t)=E[^+ f f{s,Ys,Zs,Us)ds\Tt\, F-a.s.,0<t<T. 

Step 3) We define the solution component Y of (14. ip as a progressively measurable, 
cadlag modification of 

Y{t){tu)=E[^+ f f{s,Ys,Zs,Us)ds\Tt], ujen,te[0,Tl 
where (Y, Z, U) is the limit constructed in Step 2). D 

We point out that in general under the assumptions (A1)-(A5), existence and 
uniqueness may fail to hold for bigger time horizon T or bigger Lipschitz constant 
K. See [To] for examples. However, for some special classes of generators existence 
and uniqueness may be proved for an arbitrary time horizon and for arbitrary global 
Lipschitz constants. These include generators independent of y with a delay measure 
a supported on [—7, 0] with a sufficiently small time delay 7, following Theorem 2.2 
in [in], or generators considered in [S] consisting of separate components in z and 
M, following Theorem 1 in ^. 

3 Malliavin's calculus for canonical Levy processes 

There are various ways to develop Malliavin's calculus for Levy processes. In this 
paper we adopt the approach from [23] based on a chaos decomposition in terms of 
multiple stochastic integrals with respect to the random measure M . In this setting, 
we will construct a suitable canonical space, on which a variational derivative with 

12 



respect to the pure jump part of a Levy process can be computed in a pathwise 
sense. 

In this section we give an overview of the approach of Malhavin's calculus on 
canonical Levy space according to |23] (see [23| for details). We then prove some 
technical results concerning the commutation of integration and variational differ- 
entiation, which are needed in the next section. 

We assume that the probability space {Q,J^,F) is the product of two canonical 
spaces {Qw x ^n, J^w ^ -^a^) ^w ^ '^n), and the filtration F = (J-'j)tg[o,T] the canonical 
filtration completed for P. The space {Qw, J^w, ^w) is the usual canonical space for a 
one-dimensional Brownian motion, with the space of continuous functions on [0,T], 
the cr-algebra generated by the topology of uniform convergence and Wiener mea- 
sure. The canonical representation for a pure jump Levy process (r2jv,^Af, Pat) we 
use is based on a fixed partition {Sk)k>i of M— {0} such that < i^iSk) < C)0, A; > 1. 
Accordingly, it is given by the product space ^k>ii^Nj'^N^^N) ^^ spaces of com- 
pound Poisson processes on [0, T] with intensities //(S^) and jump size distribu- 
tions supported on Sk,k > 1. Since trajectories of compound Poisson processes 
can be described by finite families {{ti,zi), ...,{tn, Zn)), where (ti,...,t„) denotes 
the jump times and {zi,...,Zn) the corresponding sizes of jumps, one can take 
^N = U„>o([0'^] X (K - {0}))", with ([0,T] X (M - {0})° representing an empty 
sequence, the a-algebra T^ = \/^yQB{{[0,T] x (M — {0}))"), and the measure P^ 
defined in such a way that for B = U„>oi?„, B^ G i3(([0, T] x (M - {0}))"), we have 

(K5.)r(rft®^)«"(5.) 



P^(E) = e-^(^^)^^ 



n! 

n=0 



Now consider the finite measure q defined on [0,T] x M by 

q{E)= / dt+ / z^u{dz)dt, EeB{[0,T]xR), 
Je{o) Je' 

where E{0) = {t G [0, T]; (t, 0) G E} and E' = E — E{0), and the random measure 
Q on [0, T]xR 

Q{E) = f dW{t) + [ zN{dt,dz), E G -B([0,T] x M). 
Je{o) J E' 

For n G N and a simple function /i„ = lEix...->cEni with pairwise disjoints sets 
Ei,...,En G i3([0,T] X M), a multiple two-parameter integral with respect to the 
random measure Q 



In{hn) = / h{{ti,Zi),...{tn,Zn))Q{dti,dZi) ■ ... ■ Q{dtn,dz„ 

J{[0,T]xR)" 

can be defined as 

Inihr.) = Q{E,)...Q{E,,). 
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The integral can be extended to the space L|. (M) of product measurable, deter- 
ministic functions h : ([0,T] x M)" — )> M satisfying 

\\h\\\2 = / \hn{{ti, zi),.. ., {tn,Zn))\'^q{dti, dzi) ■... ■q{dtn,dzn) < oo. 

The chaotic decomposition property yields that any J-'-measurable square integrable 
random variable H on the canonical space has a unique representation 



H = J2Uhn), P-a.s., (3.1) 



n=0 



with functions /i„ G Ly „(M) that are symmetric in the n pairs {ti,Zi), 1 < i < n. 
Moreover, 

oo 

E\H^]=Yn\\\hnfr2 . (3.2) 

n=0 

In this setting it is possible to study two-parameter annihilation operators (Malliavin 
derivatives) and creation operators (Skorokhod integrals). 

Definition 3.1. 1. Let ©-"^'^(M) denote the space of ¥ -measurable random vari- 
ables H G L^(]R) with the representation H = '^^=q In{hn) satisfying 



y ^nn\\\h. 



ra=l 



„||^2 < OO. 

T,q,n 



2. The Malliavin derivative DH : f] x [0,T] x R — ;> R o/ a random variable 
H G D^'^(M) is a stochastic process defined by 

oo 

Dt,zH = ^nln-iihniit,z),-), valid for q - a.e.{t, z) G [0,T] x M,P-a.s.. 

n=l 

3. Let L^'^(]R) denote the space of product measurable and ¥ -adapted processes 
G:nx[0,T]xR^R satisfying 

G{s,y) gD1'2(R), forq-a.e.{s,y) G [0,T] x M, 
^[/([o,T]xR)2 \Dt,zG{s,y)\^q{ds,dy)q{dt,dz)] < oo. 

In terms of the components of the representation ofG{s,y) = Yl'^=o-^n{gn{{s,y), ■), 
for q-a.e. {s,y) G [0,T] x M, the above conditions are equivalent to 

oo 

Y^{n + l){n + l)\\\gn\\l2 < oo, 

n=l 
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where cjn denotes the symmetrization of Qn with respect to all n + 1 pairs of 

variables. 

The space h^''^(R) is a Hilbert space endowed with the norm 



\G\\l,,=E[ [ \Gis,y)\'qids,dy)] 

JlO,T]xR 

+E[ f \Dt^,G{s,y)\'q{ds,dy)qidt,dz)]. 

J(\0.T]xR)'2 



'([0,T]xIR)2 

4- The Skorokhod integral with respect to the random measure Q of a process 
G : i7 X [0, T] X M — )■ M with the representation G{s, y) = Yl'^=o ^n{gn{{s, y), ■), 
for q — a.e.(s, y) G [0, T] x M., satisfying 



J^'^n + ly.WgJl^ <oo, 

f * T,q.n + 1 

n=0 

is defined as 

/ G{s,y)Q{ds,dy) = y2^n+i{gn), f - a.s. 

J\0,T]xR „_n 



n=0 



The following practical rules of differentiation hold. Consider a random variable 
H defined on Q^y x fijv- The derivative DffiH is with respect to the Brownian motion 
component of the Levy process, and we can apply classical Malliavin's calculus for 
Hilbert space-valued random variables. If for P^-a.e. un G Qn the random variable 
H{.,un) is differentiable in the sense of classical Malliavin's calculus, then we have 
the relation 

Dt,oH{uJw,ujN) = DtH{.,LON){cow), A-a.e.tG [0,T],P^xP^-a.s., (3.3) 

where Dt denotes the classical Malliavin derivative on the canonical Brownian space, 
see Proposition 3.5 in |21]. In order to define Df^F for -2 7^ 0, which is a derivative 
with respect to the pure jump part of the Levy process, consider the following 
increment quotient operator 

where 10^ transforms a family wat = [(t\^z\)^ (^2,-22), •••)) ^ '^n hito a new family 
uj^ = {it,z), {ti,zi), (^2,-22), •••)) ^ ^N, by adding a jump of size z at time t into the 
trajectory. According to Propositions 5.4 and 5.5 in |23], for H G L^(]R) such that 
^[/n /m-IoI I^*.^-^I f^idz)dt^ < 00 we have the relation 

Dt,,H = -^t^zH, for A m - a.e.(t, z) G [0, T] x (R - {0}), P - a.s.. (3.5) 
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The operator (13. 4p is closely related to Picard's difference operator, introduced 
in [21], which is just the numerator of (13 ■4p . It is possible to define Malliavin's 
derivative for pure jump processes in such a way that it coincides with Picard's 
difference operator, see [12]. We point out once again that we adopt the approach of 
|24) . and define multiple two-parameter integrals with respect to the random measure 
M and not with respect to A^, to obtain differentiation rules (13. 3p and (13. 5p . 

We now discuss some technical problems arising in the next section in the context 
of the main theorem of this paper. The subsequent lemmas are extensions of classical 
Malliavin differentiation rules to the setting of the canonical Levy space. 

Lemma 3.1. Assume that H G Di'2(M). Then, forO < s <T, E[H\J^s] e B^'^{R) 
and 

Dt^,E[H\Ts] =E[Dt,zH\J^s]Ht<s}, forq-a.e.{t,z) e[0,T]xR,F-a.s.. 

Proof: 

The proof is a straightforward extension of the proof of Proposition 1.2.8 from |19| . 
Details are left to the reader. D 

We next provide a proof of the commutation of Lebesgue's integration and vari- 
ational differentiability, which is commonly used. 

Lemma 3.2. Let G : Q x [0,T] x W ^ M. be a product measurable and ¥-adapted 
process, rj on [0,T] x M a finite measure, so that the conditions 

G{s, y) e ©i'2(R), for ri - a.e.{s, y) G [0, T] x M, (3.6) 

^[/([o,r]xR)2 \Dt,zG{s,y)\^ri{ds,dy)q{dt,dz)\ < oo. 

are satisfied. Then JrQrp,^G{s,y)r]{ds,dy) G ©""^'^(M) and the differentiation rule 
Dt,z G{s,y)r]{ds,dy) = Dt,zG{s,y)r]{ds,dy) 

J[0,T]xR J[0,T]xR 

holds for q-a.e. {t,z) G [0,T] xM,P-a.s.. 

Proof: 

As for rj-a.e. {s,y) G [0,T] x R the random variable G{s,y) is J-'^-measurable and 
square integrable, the chaotic decomposition property on the canonical space (13. ip 
provides a unique representation 

oo 

G{s, y) = Y, U9n{{s, y),.)), ii - a.e.(s, y) G [0, T]xR, P - a.s.. 



n=0 
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By part 3 of definition 13. 1[ the assumptions (13. 6p yield 

/ y^ nn\\\gn{{s,y),. )\\l2 r]{ds, dy) < oo. (3.7) 

J[o,r]xM^ 

For A^ G N let G^ be a measurable version of the partial sum of the first A^ + 1 
components given by 

N 

G^(s, y) = J2 U9n{{s, y),.)), 7]- a.e.{s, y) G [0, T] x R, P - a.s.. 

n=0 

We first prove that J.^rp,^^G'^{s,y)ri{ds,dy) G D-^'^(M) and the claimed differentia- 
tion rule holds. 

By canonical extension arguments, a Fubini type property holds for any single chaos 
component, and therefore 

G'^is,y)vids,dy) 

[o,r]xR 

= / X] / 9n{{s,y),{ti,Zi),...,{tn,Zn)) 

^[0,r]xR „^Q J{[0,T]xR)" 

■Q{dti,dzi)...Q{dtn, dzn)r]{ds, dy) 

= y2 / 9n{is,y),{ti,Zi),...,{tn,Zn)) 

n=0 >^([0,T]xR)" J[0,T]xR 

N 

■rj{ds,dy)Q{dh,dzi)...Q{dtn,dzn) = Y^UK) := H"" , (3.8) 

n=0 

with 



hniih, zi), ..., {tn, Zn)) = / fi'„((s, y), (ti, Zi), ..., (t„, Zn))vids, dy) 

J[0,T]xR 

for (ti, Zi), ..., (t„, z„) G ([0,T] X M)". Notice that by Cauchy-Schwarz' inequality, 
finiteness of t], the assumption (13. 7p and Fubini's theorem we obtain 



YJn?2!||/i„||^2 < oo. (3.9) 

n=l 

For any A^ G N, we have that H^ G D^'2(M), hence J^^^^^^G^{s,y)r]{ds,dy) G 
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©^'^(M), and, by linearity and definition 
Dt,.H'' = Dt,. [ G''{s,y)7]{ds,dy) 

J[0,T]xM. 

N 

= y^^ I gniis,y),{t,z),{t2,Z2,)...,{tn,Zn))'nids,dy)Q{dt2,dz2)...Qidtn,dz, 

1 JfO.TlxM 



„=1 .[0,T]xIR 

/ y2 ngniis, y), (t, z), (ta, Z2, )..., (t„, Zn))Q{dt2, dz2)...Q{dtn, dzn)r]{ds, dy) 

7[0,T]xR„^;L 



/ A.^G^ls, y)r]{ds, dy), q- a.e.{t, z) E [0, T] x 

J[0,T]xR 



'[0,T]) 

The diff'erentiation rule is proved for G^ . 
Finally, by (13. 9p we have 



E[|if^-if"?] + / E[\D,,,H''-D,,,H''\']q{dt,dz) 

M 

< V nn\\\hn\\l2 ^0, N,M^ 00. 

' ^ T.o.ti 



T,q,n 

n=N+l 



By closeability of the operator D we conclude that the unique limit H is Malliavin 
differentiable. The convergences G^ — )■ G P ® ?7-a.e. and DG^ — )■ DG W®ri® g-a.e. 
together with Lebesgue's dominated convergence theorem, justified by the first and 
third assumption in (13. 6p . give 

E[| / G''{s,y)ii{ds,dy)- [ Gis,yMds,dy)\'] 

J[0,T]xR JlO,T]xR 

+ f ml Dt,,G''is,y)r]{ds,dy)- [ Dt,,Gis,y)r]{ds,dy)\^]qidt,dz) ^ 0. 

7[0,T]xIR 7[0,T]xIR J[0,T]xR 

This implies the claimed equation. D 



We finally discuss the commutation relation of the Skorokhod stochastic integral 
with the variational derivative. 

Lemma 3.3. Assume that G : ^2 x [0, T] x M — > R is a predictable process and 
K[J,^^.^\G{s,y)\'^q{ds,dy)~\ < 00 holds. Then 

G G Li'2(M) if and only if J^^^^^^^G{s,y)Q{ds,dy) G B^'^{R). 

Moreover, if J^i^^^^^G{s,y)Q{ds,dy) G D^'2(]R) then, for q-a.e {t,z) G [0,T] x M, 

A,./ Gis,y)Q{ds,dy) = Git,z)+ [ Dt,,Gis,y)Qids,dy), F - a.s., 

J[0,T]xR J[0,T]xR 

and j,Qj,.^Dt^zG{s,y)Q{ds,dy) is a stochastic integral in ltd sense. 
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Proof: 

By square integrability of G, for q-a.e [s, y) G [0, T] x M, the chaotic decomposi- 
tion property yields the unique representation G{s,y) = Yl^=Q-^n{.9n{.{.s,y), .), Qn € 
^T qn+ii n>Q. Square integrabihty and predictabihty of G imphes that the stochas- 
tic integral L ^^ j^ G{s, y)Q{ds, dy) is well-defined in the ltd sense and the Skorokhod 
integral, which coincides under the given assumptions with the Ito integral (see The- 
orem 6.1 in j21]) can be defined by the series expansion J, rpT^^G{s,y)Q{ds,dy) = 
Yl'^=o -^n+iidn) according to Definition 3.1.4. The Skorokhod integral is Malliavin 
differentiable if and only if Yl'^=i{^ + ^){^ + I)'ll5'n|li2 < C)0, see Definition 

3.1.2. This series converges if and only if G G L^'^(R), by Definition 3.1.3. 
Following Section 6 in (24j . we can conclude that the required differentiation rule 
holds. To prove that the integral J^Qj,,^j^Dt^zG{s,y)Q{ds,dy) is well-defined in the 
Ito sense, it is sufficient to show that the integrand {u,s,y) t— )■ Dt^zG{s,y){u) is a 
predictable mapping on fi x [0, T] x M, as square integrability is already satisfied by 
G e Li'2(]R). For g-a.e. {s,y) G [0,T] x M, predictability of G implies that 

CXD OO 

G{s,y) = J2U9n{{s,y),.) = J2U9n{{s,y),.)lf,l^{.)), F - a.s., 

n=0 n=0 

and applying Definition 3.1.2 of the Malliavin derivative yields 

OO 

D,,zGis,y) = J2^In^ii9n{is,y),{t,z),.)lf,l^{{t,z),.)), 



n=0 



for q(g)q-a.e.{{t,z),{s,y)) G ([0,T] x M)^P-a.s., 

from which the required predictability of the integrand follows. As a by-product, let 
us note that [u, s, y, t, z) i— t- Dt^zG{s, y){u}) is jointly measurable. D 



4 Variational differentiability of a solution 

The main goal of this paper is to investigate Malliavin's differentiability of a solution 
of a backward stochastic differential equation with a time delayed generator. In this 
section, additionally to (A1)-(A5), we assume that 

(A6) the generator / is of the following form 

f{t,yt,zt,ut) 
■= f((^,t, / y{t + v)a{dv), 

z{t + v)a{dv) / / u{t + V , z)m{dz)a{dv] 
-T J-T Jr-{o} 
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with a product measurable function f : Q x [0,T] xMxMxM— j-M, which is 
Lipschitz continuous in the last three variables for Peg) A-a.e. {u, t) E Qx [0, T], 
more precisely the generator satisfies (A3) with the same constant K, 

(A7) the terminal value is Malliavin differentiable, i.e. ^ G D^'^(]R), and 

E[ f \DsM^q{ds,dz)] <oo, 

J[0,T]xR 

limE[ / / \Ds,z^\^m{dz)ds] = 0, 

^^^ Jo J\z\<e 



(A8) for P® A-a.e. {u!,t) E Q x [0,T], the mapping {y,z,u) h-). f(ijj,t,y,z,u) is 
continuously differentiable in {y,z,u), with uniformly bounded and contin- 
uous partial derivatives fy, fz, fu', we assume fy{u,t, ., ., .) = fz{co,t, ., ., .) = 
fu{oj, t, ., ., .) = for a; G f2, t < 0; 

(A9) for {t,y,z,u) G [0,T] x M x R x R we have f{-,t,y,z,u) G ©^'^(R) and 

E[ / / \Ds,zfi;t,0,0,0)\^dtq{ds,dy)] < oo, 

J[o,t]xrJo 

\Ds,zfiuj, t, y, z, u) - \Ds^zf{^, t, y, z, u) \ 
< L{\y-y\ + \z- z\ + \u-u\), 

(s, z) G [0, T] X R, (y, z, u) G M X R X M, {y,z,u) eRxRx M, for P ® A-a.e. 
{uj,t) eftx [0,T]. 

The assumptions (A7)-(A9) are classical when dealing with Malliavin's differen- 
tiability, see Proposition 5.3 in [13] or Theorem 3.3.1 in jTJ] in the case of BSDEs 
driven by Brownian motions. We also remark that the generator in (A6) depends 
on J_rp j]^_rQ',u{t + v,z)m{dz)a{dv), which corresponds to a standard form of de- 
pendence appearing in BSDE without delays and with jumps, see Proposition 2.6 
and Remark 2.7 in [5j. 

We can state our main theorem. 

Theorem 4.1. Assume that (A1)-(A9) hold and that time horizon T and Lipschitz 
constant K of the generator f are sufficiently small, such that for some /3 > 

1 r^ 

6 := 6{T, K, /3, a) = (ST + -)K / e-'^''a{dv) max{l, T} < 1. 

P J-T 

1. There exists a unique solution {Y,Z,U) G ^^(R) x W^iR) x M^(M) of the time 
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delayed BSDE 



+ / f(uj,r, Y{r + v)a{dv), 



t ^ J-T 

Z{r + v)a{dv), / / U{r + v, z)m{dz)a{dv)]dr 
-T J-T Jr-{q} ' 

Z{r)dW{r)- f f U{r,y)M{dr,dy), 0<t<T. (4.1) 

t Jt Jr-{o} 

2. There exists a unique solution {Y''^,Z''^,U'^^) G §2(M) x H2(M) x H^(M) of the 
time delayed BSDE 

r«'0(t) = Dsfli+ f f''\r)dr- f Z''%r)dW{r) 
Jt Jt 

(4.2) 







- [ [ U^''ir,y)Midr,dy), 
Jt Jr-{o} 


< 


s<t<T, 


with the generator 












r'% 


r) 


r\ 










= 


Dt,of 


(u,r, Y{r + v)a{dv), 












- 


/ Z{r- + v)a{dv), / / 

/-T J-T Jr- 

n 


U{r + v 

{0} 


,y)m{dy)o 


idv) 


+ 


f ( 


,r, / Y{r + v)a{dv), 












/»n 


/ Z{r + v)a{dv), / / 
J-T J-tJr 


f/(r + ^ 

-{0} 


v,y)m{dy)> 


Oi{dv 




■L 


Y''\r + v)a{dv) 










+ 


fz{0O. 


.0 

, r, / Y{r + v)a{dv), 












nC\ 


/ Z{r + v)a{dv), / / 
J-T J-T Jr. 


-{0} 


v,y)m{dy)> 


a ((if 




■L 


Z'^^{r + v)a{dv) 










+ 


f ( 


.0 

, r, / Y{r- + v)a{dv), 

pO pO p 












rtO 


/ Z{r + v)a{dv), / / 
l-T J-tJr- 


U{r + v 

-{0} 


,y)m{dy)o 


'Xdv) 




■L 


f U''\r + v,y)m{dy)a{dv). 
Jr-{o} 









(4.3) 
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3. There exists a unique solution {Y''\ Z''^,U''^) G S^{R) x B.'^{R) x H^^(M) of the 
time delayed BSDE 

Y'^'{t) = D,^,i + I f'^'{r)dr - j Z'^'{r)dW{r) 

f U'''{r,y)M{dr,dy), 0<s<t<T,z^O, (4.4) 
t Jr-{o} 

with the generator 

f(uj''\r,z f Y'''{r + v)a{dv)+ f Y{r + v)a{dv), 
Z''''{r + v)a{dv)+ I Z{r + v)a{dv), 

■T J-T 

z / U^'^{r + v^y)m{dy)a{dv) + / / U{r + v,x)m{dy)a{dv 

J-T Jr-{0} J-T Jr-{0} 

f(u,r, / Y{r + v)a{dv), 

Z{r + v)a{dv), / / U{r + v,y)m{dy)a{dv)] > /z, (4.5) 

T J-T Jr-{o} ' } 

where we set 

Y'^'{t) = Z''\t) = W'\t, y) = 0, {y, z) e (R - {0}) x M, P - a.s., t<s< T.(4.6) 

Then (Y, Z, U) G Li'2(M) x 1^^^{R) x Li'2(M) and 
(r'''^(t),Z''^(t),f/'''^(t,y))o<^,i<T,(y,2)G{K-{o})xM «s ffl version of 
(Dg^zYit), Dg^^Zit), Dg^JJit, y))o<s,t<T,(j/,z)e(R-{o})x]R- 

We recall that Dt,of{r, ., ., ., .) appearing as the first term in (14. 3 p is the Malliavin 
derivative of / with respect to u, whereas u'"^'^ appearing in (14. 5 P is defined in (13. 4p . 
Proof: 

We follow the idea of the proofs of Proposition 5.3 in [13j, or Theorem 3.3.1 in [15j . 
Let us denote by C a finite constant which may change from line to line. 
Step 1) Given /3 > 0, we prove existence of unique solutions of the equations 
(14. ip . (14. 2 p and (14. 4 p for a time horizon T and a Lipschitz constant K fulfilling 
S{T,K,f3,a) < 1. 
The existence of a unique solution {Y,Z,U) G ^^(M) x m'^{R) x e^(M) of flTT]) 



follows from Theorem 12. ![ since the assumptions (A1)-(A5) are satisfied. Under 
the additional assumptions (A7)-(A9), the time delayed BSDE (H^ and (H^ . 
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with the generators (14. 3p resp. (14. Sp . fulfill the conditions of Theorem 12.11 In par- 
ticular the corresponding generators are Lipschitz continuous with the same Lip- 
schitz constant K that the generator / possesses. It is easy to see that the gen- 
erators (14. 3 p and (14. 5 p have the same Lipschitz constant K in the sense of (A3). 
Hence, 6{T,K,^,a) < 1 holds simultaneously for all BSDEs (jH]), (g^D and (lOp 
and we conclude that for g-a.e. {s,z) G [0,T] x M there exists a unique solution 
(Y''',Z''',U''') e §2(M) X e2(M) X e^(M) of (Q or (iaD satisfying (USD- 

Step 2) Consider a sequence (F", Z", U"')n£N, constructed by Picard iteration scheme, 
which converges to {Y,Z,U). In this step we show that (y", Z", f/") G L^'^(M) x 
Li'2(M) X Li'2(M) implies (yn+i^ ^n+i^ f/n+i) ^ Li'2(M) x Li'2(M) x h^'^{R), and 
that from E[ L y, sup^gjQ^] \Ds^zy^(t)\'^(l{ds, dz)'j < oo we can as well deduce that 

^[/[o,T]xRSuPiG[o,T] \Ds,zY''+\t)\^q{ds,dz)] < oo. 



For that purpose, we study the iterations 
V^+^t) =^+ f r{r)dr 



t 
Z''+\r)dW{r)- f f U''+\r,y)M{dr,dy), < t < T, (4.7) 

where we denote 

Y''{r + v)a{dv), / Z'^ir + v)a{dv), / / W'ir + v,y)m{dy)a{dv)). 

-T J-T J-T Jr-{0} 

We first establish Malliavin's differentiability of j^ f"'{r)dr by applying Lemma 13^ 
Notice that r"(t) G Di'2(]R), for A-a.e. t G [-T,T]. Similarly to (O, we can derive 



\Y'''{r + v)\^a{dv)\dr = ¥.[ / |r"(r + i;)pdra(dv)] 

T J-T Jo 

/O j-T+v 
/ \Y''{w)\'^dwa{dv)]<TE[ sup |r"(u;)n] < oo 



oG[0,T] 

together with 



/" E[ /" /" |D,,^F"(r + t;)|2a(dt;)g(ds,d2)]dr 

< TE[ sup |D,,^y"(«;)pg(c/s,rfz)] < oo. 



'[0,T]xRu>G[0,T] 

This provides the assumptions of Lemma l3.2[ and for A-a.e. r G [0,T] we have 
J,j,Y^{r + v)a{dv) G D^'^(]R), and furthermore 

/o ^0 

Y''{r + v)a{dv)= D.^.Y^'ir + v)a{dv), F - a.s., 
-T J-T 
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for q ® A-a.e. (s, z^ r) G [0, T] x R x [0, T]. In an analogous way we derive 

Ds,J Z^{r + v)a{dv)= f D.,,,Z"{r + v)a{dv), 

Ds,z / U''{r + v,y)m{dy)a{dv)= / D,^,U''{r + v,y)m{dy)a{dv), 

J-tJr-{o} J-tJr-{o} 

holds P-a.s. for q (g) A-a.e. {s,z,r) G [0,T] x M x [0,T]. We claim that for A-a.e. 
r G [0,T] the random variable /"(r) G D^'^(]R) and for q (g) A-a.e. {s,z,r) G [0,T] x 
M X [0, T] we have 

= Dt,of{;r,j^ Y^{r + v)a{dv), 

r-O rO 



Z''{r + v)a{dv), / U'\r + v,y)m{dy)a{dv)) 



'-T J-T JIR-{0} 



/-T 

/-O 



Z"(r-Ft;)a((it;), / / U"-{r + v,y)m{dy)a{dv)\ 



-T J-TJlR-{0} 

Ds^QY'\r + v)a{dv) 

-T 

+ /.(■,r,y^ Y^{r + v)a{dv), 



Z"(r + t;)a(rft;), / / V'ir + v,y)m{dy)a{dv)) 



'T J-TJW-{Q} 

r-O 



Z''{r + v)a{dv), / / ^{r + v,y)m{dy)a{dv)) 
\ Ds,oU''ir + v,y)m{dy)a{dv), {A.l 
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and for 2; 7^ 



f(-''\r,z f Ds,,Y''{r + v)a{dv)+ I F"(r + t;)a(rft;), 

/o ^0 

Ds,,Z''{r- + v)a{dv)+ j Z'^ir + v)a{dv), 
-T ' J~T 

z / Ds^zU'^ir + v,y)m{dy)a{dv)+ I / {/'"'{r + v,x)m{dy)a{dv)) 

J-T Jr-{0} ' J-T Jr-{0} ' 

-/(■,r,y^ Y\r^v)a{dv), 

Z''{r + v)a{dv), f f W'ir + v,y)m{dy)a{dv))\ /z. (4.9) 

-T J-T Jr-{0} ' J 

The derivative (14. 8 P follows from the chain rule for the operator D^^^ as for Theorem 
2 in [22], whereas (14. 9 p follows from Proposition 5.5 in [24j provided that 

E[/ |r(r)|2dr]<oo, 

E[ [ [ [ \Ds,zrir)\^m{dz)dsdr] < 00, 
Jo Jo Jr-{o} 

hold. The finiteness of the first integral is obvious. The second integral can be shown 
to be finite by applying the Lipschitz continuity of the generator (A3) , the Lipschitz 
continuity of the derivative of the function / with respect to u and its square inte- 
grability (A9), as well as the assumption (F", Z", [/") G Li'2(M) xLi'2(M) xLi'2(M). 
Moreover, 

rT pT 



E[ [ [ \Ds,or{r)\''drds] < 00, 
Jo Jo 



-T ,„ 



and by Lemma 13.21 again we derive that for < t < T we have ^ + /^ f""'{r)dr G 
ro^'^(R) with Malliavin derivative 

Ds,.i + j DsJ'\r)dr, g - a.e.(s, z) G [0,T] x M, (4.10) 

where Dg^f^ is defined in (14. 8 p and (14. 9p . If we combine this result with Lemma 
13. H we can conclude 

F"+i(t) = E[e + I r{r)dr\J^t] e B^'\R), 0<t<T, 

and from the equation (l4.7p we derive 

f Z"+i(r)rfVr(r) gD1'2(M), 0<t<T, (4.11) 
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and 



f/"+^(r,|/)M(rfr,rfy) gD1'2(M), < t < T. (4.12) 

t JR-{0} 

Therefore Lemma [331 yields (Z"+\f/'^+i) G h'^'^iR) x h^'^{R). 

This allows us to differentiate the recursive equation (14. 7p and obtain for g-a.e.(s, z) G 

[0, T] X M 

[ Ds,M''+\r,y)M{dr,dy), s<t<T, (4.13) 
t Jr~{o} 

and 

Ds,.Y^^\t) = D,^,Z^+\t) = Ds,,U^+\t,y) = 0, t<s,ye(R- {0}). (4.14) 
Note that the time delayed BSDE JKW) with generator (gS]) or (gH) fulfills the as- 



sumptions of Theorem 12.11 with zero corresponding Lipschitz constant. We conclude 

that for g-a.e.(s, z) G [0, T] xR there exists a unique solution (D,,^F"+^ D,,^Z"+\ L),,^t/"+i) 

§2 X H^ X H^ of IKT^ satisfying IKU^ . By applying Lemma I2II1 with f = and 

f = f, together with the estimate (12.71) . with 6 = 6(T,K, f3,a) < 1 we derive the 

inequality 

\\p) ■i/'n+l||2 illn 7n+l||2 iMt-) rr"+l||^ 

< 9e^''E[\Ds,M^]+6{\\D,^,Y%, + \\D,^,Z%, + \\Ds,zU%,J. (4.15) 

This in turn yields E[ L^, jgSup^g^o.T] \Ds,zy^^^i't)\'^Qids, dz)~\ < oo, and in partic- 
ular, F"+i G Li'2( 



Step 3) We establish the integrability of the solution Y^'^(t), Z'^'^{t),U^'^{t,y) with 
respect to the product measure q on ([0,T] x M)^. 

Take {s,z) G [0,T] x M. Consider the unique solution {Y'^', Z''^,U''^) G S^{R) x 
H2(M) X H^(M) of the equation (g^D or (i3D- Lemma EH with f = and / = / 
together with the estimates (12. 7p and (12.81) yield the inequality 

II ii2 II ii2 II ii2 

r ' —I— / ' —I— / / ' 

Ir ll§2 "^ 11^ ||e2 "^ ir Ile2 



\ »C R^\^\l^s,z<i\ J T L'Vll-' ||g2 T 11^ ||jj2 T ||U ||jj2 ^5 

SO that under 6 = 6{T, K, /3, a) < 1 we obtain for g-a.e. (s, z) G [0, T] x M 

r^lli + Il^^1le2 + Pnk ^ Cn\DsM'l (4.16) 
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and we arrive at 



E[ f \Y'''{t)W{dt,dy)q{ds,dz)] < oo, 

i([0,T]xR)2 

E[ f \Z'^'{t)W{dt,dy)q{ds,dz)] < oo, 

i([0,T]xR)2 

E[ f \U'^'{t,y)\'^q{dt,dy)q{ds,dz)] < oo. 

i([0,T]xR)2 

Step 4) We show convergence of (F", Z", t/")„eN in 1L^'^(R) x Li'2(M) x L^'2(M). 
From Theorem 12.11 we aheady know that (F", Z", f/'^)„gN converges in S^(]R) x 
EI^(M) X EI^(]R). We have to prove that the corresponding Malhavin derivatives 
converge. The convergence 



J[0,T]xR 

+ \\Z''' - D,,,Z"+i||^, + \\U''^ - /},,,[/"+! 11^, )q{ds, dz) = 0, 



for z = can be proved in the similar way as in the case of a BSDE without delay 
driven by a Brownian motion, see for example Theorem 3.3.1 in [15]. We only prove 
the convergence for 2; 7^ 0. 

Lemma EB applied to the time delayed BSDE (113]), and <KW^ with (HHD, yield 
the inequality 

\\l — 'LJs,zJ- ||§2 + 11^ — -L^s,z^ Ile2 "*" 11'^ ~ J^s,z'^ ||jj2 

<{8T+^)E[J e^^\r^^ir)-D,,J-ir)\'dr], (4.17) 

for g-a.e. {s,z) G [0,T] x R. 

First, by the Lipschitz continuity condition (A3) for the generator / and the Lips- 
chitz continuity condition (A9) for the derivative of the function / with respect to 
u we obtain for A ® m ® A-a.e. (s, z, r) G [0, T] x (M — {0}) x [0, T] the following two 
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estimates 



ni^\\2 



|r'^(r)-D,,,r(r)| 



< C(f \Y'^'ir + v)\^aidv)+ f \Z'^'ir + v)\^a{dv) 
+ / \U''''{r + v,y)\'^m{dy)a{dv) 

^0 

|D,,,F"(r + v)\'a{dv) + / |/^.,.Z"(r + v)\'a{dv) 

T J ~T 

. 

/ \D,^,U'\r + v,y)\''m{dy)a{dv) 

-T JR-{0} 

^0 

+ / |F"(r + t;) - r(r + v)\^a{dv) + / |Z"(r + v) - Z{r + t;)|^a(rft;) 

-T J-T 

r-0 



|?7"(r + t;, y) - f/(r + v, y)\^m{dy)a{dv)) , (4.18) 

and for any A > 

1 / /"^ 

+ / \Z'''{r + v)- Ds,,Z''{r + v)\'^a{dv) 
r 

2. 



+ / / \U''\r + v,y)-D,,,U''{r + v,y)\'midyHdv] 

J-T Jr-{0} 

+ (1 + A) (2 + j)k(^ I \Y{r- + v)- Y^{r + v)\^a{dv) 

/■O 

+ / \Z{r\v)-Z'\r\v)'^a{dv) 

+ [ f \U{r + v,y)-U''{r + v,y)\'^m{dy)a{dv))/z'^. (4.19) 



-T jR-{0} 

Note that 

/ \\\J: ^s,z-i: ||g2 

^[0,T]x(E-{0}) 

+ ||Z^'^ - I^s,.^"+'||e2 + \\U''' - Ds,M''^%, )qid.s, dz) 



|2 



lim [ [ (||r"'"-D,,r"+H 

+ ||Z^'^ - D,,,Z^^%, + ||f/^'^ - D,,,f/"+i|£, )m(rf^)rfs. (4.20) 
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We prove that this convergence is uniform in n. 

Choose 6 > sufficiently snialL By assumption (A7) we can find e such that 



E[ / / \DsM^m{dz)ds] < e, 

Jo J\z\<e 



and 

m{dz) < e. 



z\<t 



Take arbitrary < ei < 62 < e. By applying the inequality fl4.17p . the estimate 
fl4.18p and by similar calculations as in fl2.7p we can derive 



^0 JeT<\z\<eo 



^n+l||2 

...52 
'0 Je^<\z\<e2 

-T r rT 



S 



< C f f E[ f e^'\f'''{r)-Ds,,rir)\'^dr]midz)ds 

Jo JeT<\z\<eo. Js 



< C\ I I {\\Y-f^, + \\Z-f^, + \\U-f^,)m{dz)ds 



Jei<\z\<e2 
T 



+ [ f {\\D,,,Y%, + \\Ds,.Z%, + \\D,,M%, )m{dz)ds 

Jo Jei<\z\<e2 

^0 Jei<\z\<e2 " -" 

(4.21) 

To estimate the first term on the right hand side of (I4.2ip . notice that the inequality 
(HTTB]) yields 



^0 Jti<\z\<e2 

<CE[f f \D,^,^\'^m{dz)ds] <Ce. (4.22) 

Jo Je-,<\z\<e9 
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Recalling 6 = 6{T, K, /3, a) < 1 and applying the inequality (14.151) we estimate the 
second term in fl4.2ip by 



[ [ i\\^s,zY%, + \\Ds,.Z%, + \\Ds,M%, )m{dz)ds 

Jo Jei<\z\<e2 "" 

< 9e^^E[ f f \DsM'^m{dz)ds] 

Jo Jei<\z\<e2 

+S [ f (||Z^,,.r"-i||g, + ||Z^,,,Z"-i|£, + \\Ds,.U'''-%, )midz)d. 

Jo Je,<\z\<€o, 



< 



l-S 



+6^ [ [ {\\Ds,zY% + \\Ds,zZ%2 + \\Ds,zUm, )m{dz)ds. (4.23) 

Jo Jei<\z\<e2 "" 

The estimate of the third term follows from the contraction inequality (12.81) 



Jo Jei<\z\<e2 

< 5"T(||rO - r 11^, + 11^0 - Z\\l, + ||f/° - U\\l, )e. (4.24) 



Choosing Y^ = Z^ = U^ = {) and combining (I4.22p - (l4.24p gives the uniform conver- 
gence of (i2DD- 

Next, by applying the inequality (I4.17p . the estimate (I4.19p and similar calculations 
as in (12. 7p and (12. 8p we can derive 



f [ {\\Y'^' - Ds,zY'' 

Jo J\z\>e 



An+l||2 
.Is2 
'O J\z\>e 

+ ||Z^'^ - Ds,,Z^+%2 + \\U''' - Ds,,U"+%, )m{dz)ds 



{ST+h [ [ E[ [ e^'-ir'^r) ~ D,^J^{r)\^dr]m{dz)ds 

P Jo J\z\>e Js 

*- ^ Jo J\z\>e 

+ \\Z''' - Ds,zZ%2 + \\U''' - Ds,zU^\\l2 )m{dz)ds 

1 1 ini 1 1 -"^m 

+ (1 + ^) (2 + \) (r- - y\\l. + \\Z- - Z\\l2 + \\U- - U\\lj I u{dz)]^ 

and we choose A sufficiently large such that 5 := 5(l + j) < 1. 

Due to the convergence of (F", Z", f/")„gN, for a sufficiently small £ > we can find 

A^ sufficiently large such that for all n > A^ 

(1 + A) (2 + \) {\\Y- - y 11^, + ||Z" - Z\\l + \\U-- UCJ I uid 



z] < e. 

\z\>e 
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We derive the recursion for n > N 



I I (\\\^s,z _ p) -v^n+l||2 

/ / Ur -^s,2:-f ||g2 

+ ||Z^'^ - Ds,,Z^+'\\l, + \\U-'^' - Z^,,,^"+i|£, )midz)ds 

II 1 1 ini 1 1 II inijj^ 

< ~5{ [ [ {\\Y-'^^ - D,,,Y% 

Jo J\z\>e 
Jo J\z\>e 

+ \\Z''' - Ds,zZ^\\l, + \\U''' - D,,,U''\\l,Jm{dz)ds + 
and finally we conclude that 

hm / / {\\Y^'' - Ds,.Y-+'\\l 

II 1 1 io. 1 1 II ini;,^^ 

The equation 



Se 



1-5' 



lim / (||F^'^-D,,,r"+i||g, 

"■^°° J[0,T]x(]R-{0}) 

+ \\Z''' - Ds,,Z^+%, + ||[/^'^ - Ds,,U"+'\\l, )q{ds,dz) = 

II 1 1 io. 1 1 II inijj^ 

now follows by interchanging the limits in n and e in (I4.20p . 

Step 4) Since the space L^'^(M) is a Hilbert space and the Malliavin derivative is 
a closed operator, see Theorem 12.6 in [T2], the claim that {Y,Z,U) G L^'^(M) x 
Li'2(M) X Li'2(M) and {Y-'^'{t), Z'^%t),U-'^%t,y))o<s,t<T,iy,z)m-{om is a version of 
the derivative {Ds,zy(t), D^^^Zit), Ds,zU{t,y))o<s,t<T,iy,z)eiR-{o})R follows, and fin- 
ishes the proof. D 

The following Corollary shows that the interpretation of the solution component 
{Z, U) in terms of the Malliavin trace of Y still holds for BSDE with time delayed 
generators. 



Corollary 4.1. Under the assumptions of Theorem\4-l\ we have 



((A,o>^)^(i))o<t<T *^ « version of {Z{t))^^^^^, 
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where (■)'^ denotes the predictable projection of a process. 

Proof: 

The solution of ( 14. ip satisfies 

Yis) = r(o) 

-jf{rj Y{r + v)a{dv), 

Z{r + v)a{dv), / / U{r + v,y)7n(dy)dv]dr 

-T J-T Jr-{o} ' 

+ Z{r)dW{r)+ / U{r,y)M{dr,dy), 0<s<T. (4.25) 
Jo Jo Jr-{o} 

By differentiating fl4.25p we obtain according to Lemma I3l3] for g-a.e. {u, z) G [0, T] x 



D^flY{s) = Z{u) 



r D^^ofir)dr+ f D^^,Z{r)dW{r) 

u J u 

+ / / Du,oU{r, y)M{dr, dy), 0<u<s<T, 
'u Jr-{o} 



and for z 7^ 

D„,,y(s) = U{u,z)~ f D^^J{r)dr+ f D^^,Z{r)dW{r) 

J u J u 



+ / / Du^,U{r, y)M{dr, dy), < u < s < T, 

Ju Jr-{o} 

where the derivative operators D^^z are defined by (14. 3 1) and (14. 5p . Since the mappings 
s i-> J Du^zf{r)dr,s I— )■ J Du,zZ{r)dW{r) are P-a.s. continuous and the mapping 

"^ ^ /m lR-io\ ^^^^^(^^y)^(^^^ ^y) ^^ P-a-s. cadlag (see Theorems 4.2.12 and 4.2.14 
in |i4j), taking the hmit s ^u yields 

D.^^qY{u) = Z{u), for A - a.e.u G [0, T], P - a.s., 

Du,zy{u) = U{u,z) ioT X^m-a.e.{u,z) e[0,T]x {R- {0}),F-a.s.. 

As y G §^(]R) has P — a.s. cadlag F-adapted trajectories, for < m < T we have the 
representation 

oo oo 

Y{u) = Y,ln{gn{{u,0),.) = ^/„(^„((n,0),. )!§;„](.)), g„. G 4,,„+i,n > 0, 

n=0 n=0 

with cadlag mappings u h-> gn{{u, 0), .). By Definition 3.1.2 of the Malliavin deriva- 
tive we arrive at 

oo 

D^,zY{u) = J2^Ir.-i{9n{{u,0),{u,z),.)lll^{{u,z),.)), ioT q - a.e.{u, z) E [0,T] X B 

n=0 

32 



For 5{o} X m-a.e. z G M, we conclude that the mapping {u,u) i— t- Du,zy{u){u}) is 
F-adapted and measurable and has a progressively measurable (optional) modifi- 
cation. Moreover, notice that the optional process u i— )■ Du^zY{u) and its unique 
predictable projection u H- (D„,^y)^('u) are modifications of each other, see Theo- 
rem 5.5 in f]3]. Finally, we remark that there exists aV x i3(M) measurable version 
of {uj,u,z) h-> {Du^zY)'^{u){u)), see Lemma 2.2 in \l\. This completes the proof. D 
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